Energy spectra of cosmic rays accelerated at ultrarelativistic shock waves 
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Energy spectra of particles accelerated by the first-order 
Fermi mechanism are investigated at ultrarelativistic shock 
waves, outside the range of Lorentz factors considered pre- 
viously. For particle transport near the shock a numerical 
method involving small amplitude pitch-angle scattering is 
applied for flows with Lorentz factors 7 from 3 to 243. For 
large 7 shocks a convergence of derived energy spectral indices 
up to the value tJoo ~ 2.2 is observed for all considered turbu- 
lence amplitudes and magnetic field configurations. Recently 
the same index was derived for 7-ray bursts by Waxman [As- 
trophys. J. Lett. 485, L5 (1997)]. 

In currently favored gamma-ray burst (GRB) models 
optically thin emitting regions move relativistically, with 
Lorentz factors of order of a few hundreds (cf. a re- 
view [1]). The power-law form of the spectrum often 
observed at high photon energies suggests the existence 
of nonthermal population of energetic particles. It was 
also proposed that GRB sources may produce cosmic ray 
particles with extremely high energies [2]. Thus mod- 
eling of burst sources requires a discussion of particle 
acceleration processes, possibly the Fermi acceleration at 
ultrarelativistic shock waves. 

The work of Kirk and Schneider [3] opened the prob- 
lem of cosmic ray acceleration at relativistic shock waves 
for quantitative consideration. Substantial progress since 
that time (e.g. Heavens and Drury [4], Kirk and Heav- 
ens [5], Begelman and Kirk [6], Ostrowski [7], Bednarz 
and Ostrowski [8]; for a review see Ostrowski [9] and 
Kirk [10]) clarified a number of issues related to shock 
waves with velocities reaching 0.98c or the Lorentz fac- 
tor 7 ~ 5 , but - to our knowledge - no one has attempted 
to discuss particle acceleration at shocks moving with ul- 
trarelativistic velocities characterized with large factors 
7 >> 1. 

The main difficulty in modelling an acceleration pro- 
cess at shocks with large 7 is the fact that involved parti- 
cle distributions are extremely anisotropic in shock, with 
the particle angular distribution opening angles ~ 7"^ 
in the upstream plasma rest frame. When transmitted 
downstream the shock particles have a limited chance 
to be scattered so efficiently to reach the shock again, 
but the energy gain of any such "successful" particle can 
be comparable to its original energy. As pointed out by 
Bednarz and Ostrowski [8] any realistic model of parti- 
cle scattering at magnetohydrodynamic turbulence close 
to the relativistic shock cannot involve large-angle point- 
like scattering. The choice is either to integrate exactly 
particle equations of motion in some "realistic" structure 



of the perturbed magnetic field, or to use a small-angle 
scattering model for particle momentum. With the an- 
gular scattering amplitude Afl << 7^^ and the mean 
scattering time At not too short [At > Tg (A51)^ , where 
Tg is the particle gyration period] , the last model repro- 
duces the pitch-angle diffusion process at small amplitude 
waves. We prefer that approach to the exact integration 
of equations of motion of a particle because of its rela- 
tive simplicity. It is also suggested that it can be reason- 
ably used for modeling particle trajectories in turbulent 
fields with large amplitude, if small At [< Tg (Afi)^] is 
involved. Below, a hybrid method involving the small 
amplitude pitch-angle scattering is applied for a particle 
transport near the shock for flows with Lorentz factors 7 
from 3 to 243. 

/. Numerical simulations. - In the present consider- 
ations we model a process of cosmic ray particle accel- 
eration by applying the following Monte Carlo approach 
(cL e.g. [8]). Energetic seed particles are injected at 
the shock and each particle trajectory is followed using 
numerical computations until it escapes through the free 
escape boundary placed far downstream from the shock 
or it reaches the energy larger than the upper energy 
limit at Emax = IO^^-Eq (-E'O - mitial energy). Simu- 
lations are continued until one obtains the power-law 
spectrum in the full range {Eq, Emax) (cf- [7,8]). All 
computations are performed in the respective - upstream 
or downstream - plasma rest frame. Each time when the 
particle crosses the shock its momentum is Lorentz trans- 
formed to the respective plasma rest frame and, in the 
shock normal rest frame (cf. [6]; henceforth "shock rest 
frame" ) the respective contribution is added to the given 
momentum bin in the particle spectrum. In the sim- 
ulations we use a simple trajectory splitting technique. 
All particles are injected into simulations with the same 
initial weight factors 1.0 . When some particles escape 
through the boundary, we replace them with ones arising 
from splitting the remaining high-weight particles in a 
way to preserve its phase space coordinates but ascribing 
only a half or a smaller respective part of the original 
particle weight to each of the resulting particles. For 
any particle crossing the shock a factor is added to the 
simulated spectrum in the shock rest frame equal to the 
particle weight divided by its velocity component normal 
to the shock. 

Efficient particle scattering with a very small AfJ re- 
quires derivation of a large number of scattering acts and 
the respective numerical code becomes extremely time- 
consuming. In order to overcome this difficulty in the 
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present simulations we propose a hybrid approach in- 
volving very small Afii (<< 7~^) close to the shock, 
where the scattering details play a role, and much larger 
scattering amplitude Af72 = 9° to describe particle dif- 
fusion further away from the shock. The respective scal- 
ing of the scattering time At is performed in both cases 
(Af2^/Ati = Ar22/Ai2) to mimic the same turbulence 
amplitudes measured by the values of the cross-field diffu- 
sion coefficient, k±, and the parallel diffusion coefficient. 

For a few instances we checked the validity of this ap- 
proach by reproducing the results for the small AOi ev- 
erywhere. In the present simulations we assume the same 
scattering conditions upstream and downstream from the 
shock (the same k± and k\\ in the units of r^c, where rgC 
is the particle gyration radius in the unperturbed back- 
ground magnetic field) , preserving particle energy at each 
scattering in the plasma rest frame. In the simulations 
we considered a few configurations of the upstream mag- 
netic field, with inclinations with respect to the shock 
normal being V = 0°, 10°,20°, 30°,60° and 90°. The 
first case represents the parallel shock, the second is for 
the oblique shock - subluminal (i.e. with the shock ve- 
locity projection at the magnetic field with slower than 
light velocity) at 7 = 3 and a superluminal one at larger 
7, and the larger ip are for superluminal perpendicular 
shocks for all velocities. The downstream magnetic field 
is derived for the relativistic shock with the compression 
R obtained with the formulas of Heavens and Drury [4] 
for a cold (e, p) plasma - i? « 3.6 for our smallest value 
of 7 = 3 and tends to ii = 3 for 7 >> 1, as measured in 
the shock rest frame. 

II. Results. - Particle spectral indices were derived for 
different mean magnetic field configurations, measured 
by the magnetic field inclination ip with respect to the 
shock normal in the upstream plasma rest frame, and for 
different amounts of turbulence measured by 

In successive panels in Fig. 1 the energy spectral in- 
dices, (J, for varying ip and k±/k^^ are presented. For a 
parallel shock {tjj = 0°) the amoimt of scattering does 
not influence the spectral index and for growing 7 it ap- 
proaches (Too — 2.2. One may note that essentially the 
same limiting value was anticipated for the large-7 par- 
allel shocks by Heavens and Drury [4]. The results for 
ip = 10° are for superluminal shocks if 7 > 5.75. 

In this case, when we go from the "slow" 7 = 3 shocks 
to higher 7 ones, at first the spectrum inclination in- 
creases {a grows), but at large; 7 the spectrum flattens 
to approach the asymptotic value close to 2.2. The spec- 
trum steepening phase is more pronounced for small am- 
plitude perturbations (small k±/k\\), but even at very 
low turbulence levels the final range of the spectrum flat- 
tening is observed. For larger ifj the situation does not 
change considerably, but the phase of spectrum steepen- 
ing is wider, involving larger values of a and starting at 
smaller velocities, below the lower limit of our consider- 



ations (there may be no such range involving the small 
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FIG. 1. The simulated spectral indices a for particles ac- 
celerated at shocks with different Lorentz factors 7. Results 
for a given kx /k|| arc joined with lines; the respective value of 
logj^Q k±/k\\ is marked by the point shape (see upper panel). 
The results for different magnetic field inclinations ip are given 
in the successive panels: (a) ■0 = 0°, (b) V = 10°, and (c) 
V' = 90°. 

steepening phase if the required velocity is below the 

sound velocity). The spectral indices for different 
magnetic field inclinations, but for the same value of 
logiQ(K;j_/«;||) = —3.44, are presented at Fig. 2. 

The large spectral indices occurring in the steepening 
phase are usually interpreted as a spectrum cutoff. In 
this case the main factor increasing the particle energy 
density is a nonadiabatic compression in the shock [6]. 
The particle angular distributions f (/u) in the considered 
shocks can be extremely anisotropic when considered in 
the upstream plasma rest frame. However, when pre- 
sented in the shock rest frame the distribution is always 
"mildly" anisotropic. This feature is illustrated at Fig. 3 
for 7 equals 3 or 27 (note that in Figs. 3-6 the area be- 
low each curve is normalized to 100). In simulations we 
observed an interesting phenomenon accompanying pre- 
viously discussed spectrum convergence to the limiting 
inclination: spectra close to the limit exhibit similar an- 
gular distributions at the shock as measured in the shock 
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FIG. 2. The simulated spectral indices a for particles ac- 
celerated at shocks with different Lorentz factors 7. Results 
for a given upstream magnetic field inclination ijj are joined 
with dashed lines; the respective value of tp is given near each 
curve. The value A = \og-^Q{K± / ku) is given in the figure. 
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FIG. 3. The simulated particle angular distributions in 
the shock in different coordinate frames: UP: the upstream 
plasma rest frame, DOWN: the downstream plasma rest 
frame, and SHOCK: the shock rest frame. The results are 
presented for parallel shocks with the Lorentz factors 27 and 
3 given near the respective curves. In the upper panel the left 
axis is for 7 = 27 and the right one for 7 = 3. 

rest frame (Fig. 4). Again, this feature is independent of 



the background conditions, and the difference between 
the actual angular distribution and the limiting one re- 
flects the difference between the spectral index a and 
CToo (cf. Fig. 5). For parallel shocks with 7 > 9, where 
the spectral index is essentially constant a = Uoo, this 
distribution is independent of the value of 7 and the per- 
turbation amplitude (Fig. 6). 

///. Discussion. - For large 7 shocks we observe the 
convergence of the derived energy spectral indices to the 
value (Too ~ 2.2, independently of background conditions. 
This unexpected result providing a strong constraint for 
the acceleration process in large 7 shocks requires a more 
detailed analysis. Particularly, the acceleration process 
in the superluminal shocks has to be clarified in the pres- 
ence of the small amplitude turbulence. 

The inspection of particle trajectories reveals a simple 
picture of acceleration. Cosmic ray particles are wander- 
ing in the downstream region with the shock wave moving 
away with the mildly relativistic velocity « c/3. Some 
of these particles succeed to reach the shock, but then 
they remain in the upstream region for a very short time 
- being very close to the shock - due to large shock ve- 
locity « c. This scenario is essentially equivalent to the 
picture involving particles reflecting in a nonelastic way 
from the receding wall. 

For large 7 shocks any particle crossing the shock up- 
stream has a momentum vector nearly parallel to the 
shock normal (cf. Ostrowski [7]); e.g., for 7 = 243 
the momentum inclination must be smaller than Omax ~ 
0.24°. If the scattering or the movement along the curved 
trajectory increases this inclination above the mentioned 
limiting value the particle tends to recross the shock 
downstream. One should note that even a tiny - compa- 
rable to Omax " angular deviation in the upstream plasma 
{A9u) can lead to large angular deviation for 7 » 1 
as observed in the downstream rest frame. The phe- 
nomenon of decreasing a to a^c at constant A and for 
growing 7 results from slower diminishing of the part of 
A6u caused by scattering in comparison to A0u arising 
due to trajectory curvature in the uniform field compo- 
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FIG. 4. Examples of the shock rest frame particle angular 
distributions for different cases with a close to (Too- 
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FIG. 5. The shock rest frame particle angular distributions 
for 7 = 27 and ■0 = 30°. Curves are presented for increasing 
A = logjQ K^l K\\ and a approaching (Too. The last curve is the 
same as the curve (b) at Fig. 4. 

nent. In this way the magnetic field structure defined 
by ■(/; becomes unimportant, at least for 7 — > oo (one 
should also note that the downstream field inclination 
approaches 90° if 1/; 7^ and 7 3> 1). As a result par- 
ticles crossing the shock downstream are scattered in a 
wide angular range with respect to the shock normal, 
providing some particles with the trajectory phase pa- 
rameter allowing for recrossing the shock upstream even 
for the perpendicular magnetic field configuration. Our 
interesting finding not fully explained with such simple 
arguments is of the belief that the resulting spectral in- 
dex is the same for oblique and parallel shocks. We also 
observe that when approaching the limiting value of the 
spectral index the mean particle energy gain < l^E/ E > 
in the cycle "upstream-downstream-upstream" reaches a 
value close (slightly above) to 1.0, much smaller than a 
factor 7^ expected for a model involving a large angle 
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FIG. 6. The shock rest frame particle angular distributions 
for parallel shocks with 7 = 3, 9, 27, 81, 243. A visible 
deviation of the distribution for 7 = 3 (full line) results from 
the slightly larger compression occurring in such a shock. 



can be roughly estimated as a fraction of the gyration 
time in this region. 

IV. Final remarks. - The presented results are to 
be applied in models of GRB sources involving ultrarel- 
ativistic shock waves. One should note that the mean 
downstream plasma proton energies can reach there sev- 
eral tens of GeV (cf. Paczynski and Xu [11]) and the 
lower limit of the considered cosmic ray energies has to 
be larger than this scale. For shocks propagating in (e~, 
e"*") plasma the involved thermal energies are lower, ~ 7 
MeV. These estimates provide the respective lower limits 
for the accelerated cosmic ray particles. For the physical 
conditions considered in GRB sources the acceleration 
process can provide particles with much larger energies, 
limited only by the condition that the energy loss pro- 
cesses (radiative, or due to escape) are ineffective in the 
downstream gyroperiod time scale. We note a striking co- 
incidence of our limiting spectral index with the value 2.3 
±0.1 derived for energetic electrons from gamma-burst 
afterglow observations [12]. 

Our derivations are limited to the test particle ap- 
proach. However, as the obtained spectra are charac- 
terized with a > 2.0, any nonlinear back reaction effects 
are not expected to affect the acceleration process within 
the spectrum high energy tail with a « (Too- 
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pointlike scattering. Thus the particle acceleration time 
scale, as measured in the downstream plasma rest frame, 
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